Quantum Poincare Recurrences 
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We show that quantum effects modify the decay rate of Poincare recurrences P(t) in classical 
chaotic systems with hierarchical structure of phase space. The exponent p of the algebraic decay 
P{t) oc l/t p is shown to have the universal value p = 1 due to tunneling and localization effects. 
Experimental evidence of such decay should be observable in mesoscopic systems and cold atoms. 

PACS numbers: 05.45.+b, 03.65.Sq 



The general structure of classical phase space in chaotic 
Hamiltonian systems displays a hierarchical mixture of 
integrable and chaotic components down to smaller and 
smaller scales [Q. This complicated structure leads, in 
particular, to an anomalous power law decay of Poincare 
recurrences P(t) and correlations C(t) inside the chaotic 
components [§1). Physically, such slow decay appears 
due to a decrease, down to zero, of the diffusion rate for 
a trajectory when it approaches the chaos border deter- 
mined by some critical invariant curve |^-||. Typically, 
P(t) ~ C(t)/t ~ t-P with p « 1.5 [§. As a result the 
integrated correlation function, which determines the dif- 
fusion rate (D ~ J Cdt), can diverge thus leading to a 
super-diffusive propagation 0. Such effects are impor- 
tant for electron dynamics in super-lattices where usually 
the phase space has a mixed structure || . 

The above anomalous properties had been studied in 
great detail for classical systems [||-||. However, the 
question how they are affected by quantum dynamics 
was not addressed up to now. This problem becomes 
more and more important not only due to its fundamen- 
tal nature but also in the light of recent experiments with 
mesoscopic systems. Indeed, different types of ballistic 
quantum dots can now be studied in laboratory experi- 
ments |!| and the phase space in such systems generally 
has a mixed structure. Since, the probability to stay in 
a given region is directly related with P(t) and C(t), its 
slow decay can significantly affect conductance proper- 
ties. In particular it has been proposed that such decay 
should lead to fractal conductance fluctuations 

@, the 

experimental observation of which has been reported re- 
cently [jy]. According to fl(],[ll| the fractal exponent a 
for conductance fluctuations is directly related to the ex- 
ponent p as (7 — 2 — p/2. 

A different type of systems in which such effects should 
be observable experimentally is given by cold atoms in 
external laser fields where the Kicked Rotator model of 
quantum chaos has been built experimentally E , ff3| . 
Possibilities of experimental investigation of slow proba- 
bility decay in such systems has been discussed recently 



s 

The experimental studies of slow power correlation de- 
cay in the regime of quantum chaos are also important 
from the fundamental point of view, since here the typical 
scale of correlation decay is much larger than the Ehren- 
fest time scale tE ~ In l/7i on which the minimal coherent 
wave packet spreads over the avaible phase space. To the 
best of our knowledge the comparison of classical and 
quantum correlations in such regime has not been inves- 
tigated so far. Only recently such comparison has been 
made in the regime of hard chaos with exponential cor- 
relation decay [|l6| . In this paper we address directly the 
comparison between classical and quantum Poincare re- 
currences (QPR), related to the correlations decay, in the 
regime with mixed phase space. Our results demonstrate 
a new universal law for QPR related to localization and 
tunneling effects. 

To investigate the QPR we use the model of kicked 
rotator with absorbing boundary conditions studied in 
|p^ , ^6| . The evolution operator over the period T of the 
perturbation is given by 



-iTfi 1 4 -ik cos6 - 



(1) 



where V is a projection operator over quantum states 
n in the interval (— N/2, N/2). Here, we put K = 1 so 
that the commutator is [n, 0] — —i and the classical limit 
corresponds to k — * oo, T — > while the classical chaos 
parameter K = kT remains constant. In the classical 
limit the dynamics is described by the Chirikov standard 
map: 

T -]re = e + T - 

2 J ' 2 
in which orbits are absorbed outside the interval —N/2 < 
n < N/2. In order to study the classical and quantum 
survival probability P(t) we fixed the ratio N/k — 4 and 
take the classical chaos parameter K = 2.5, so that the 
classical phase space has a hierarchical structure of inte- 
grable islands and chaotic components. 

A typical example of classical and quantum survival 
probability decay is shown in Fig. 1 . The classical proba- 
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bility P(t) decays with a power law with exponent p w 2 
in a range of six orders of magnitude. In this case the 
exponent is slightly different from the typical value 1.5. 
Indeed as discussed in ®^ the exponent can vary from 
system to system (and even oscillates with kit) depend- 
ing on the local structure of phase space in the vicinity 
of the critical boundary invariant curve whose rotation 
number can play an important role. 

The quantum survival probability P q (t) is plotted for 
different values of N ( effective h is proportional to l/N) 
in Fig.l. In agreement with the correspondence princi- 
ple, the quantum probability follows the classical value 
during a rather long time scale which grows with N. For 
longer times the quantum probability P q (t) is found to 
approximately follow the quantum decay law 

P q {t) « C/t (3) 

where C is some constant. This latter decay continues 
during a rather long interval of time (4 orders of magni- 
tudes in time for the case N — 3 6 ). Of course (g) repre- 
sents an intermediate asymptotic behavior since for t — ► 
oo the decay will be exponential P q (t) ~ exp (— Y m i n t) 
where Y m i n is determined by the minimal imaginary 
part of the eigenvalue A of the evolution operator U 
(J7Va = e- iX ip x ). 

In order to better understand the origin of quantum be- 
havior (|^) we studied the probability decay in a simpler 
case where the classical dynamics is completely chaotic 
(K = 7) or quasi-integrable (K = 0.5). The quantum 
and classical probability decays are shown in Fig. 2 for 
k = 5 and absorption for n < and n > N — 500. Ini- 
tially the probability is concentrated at n — 0. For K = 7 
the classical probability decays, asymptotically, exponen- 
tially with time P(t) = 0.11 exp (—"ft). The value of 
7 can be found from the solution of the Fokker-Planck 
equation with absorbing boundary conditions which gives 
7 w Dn 2 /2n 2 , where D — (3k 2 / 2 is the classical diffusion 
rate (see eq.(4) in |15[). The value of (3 depends on the 
classical chaos parameter and for K — 7 is f3 ~ 2.8 [ PHI . 
Therefore the expected theoretical value is 7 = 6.9 • 10~ 4 
which is close to the numerical value 7 = 6.4 • 10~ 4 . The 
asymptotic exponential decay starts after the diffusive 
time to ~ 1/7- The quantum probability in this case 
(Fig. 2) follows the classical one during some interval of 
time after which it decays according to eq. (||) . 

The quantum decay ([5]) can be understood in the fol- 
lowing way. The quantum eigenstates are localized with 
localization length £ D/2 [Tq] . Therefore the absorp- 
tion time for a state at a distance n from the absorbing 
boundary is t ~ exp (2n/£). Since n is proportional to the 
total measure fi ~ n oc £lnt, it follows that the survival 
probability is P q (t) = dfi/dt ~ £/t in agreement with (|^). 
The same estimate can be obtained by expanding an ini- 
tial state over eigenstates of the evolution operator (|l|) 
with probabilities |c„| 2 ~ e ~ 2n / 1 /£ and ionization rates 



T„ ~ |c„| 2 ~ e 2n / £ /D. Then the survival probability 
P q (t) ~ J °° |c„| 2 e- T ^dn ~ D/t. Our data at K = 7 for 
different values of k (k — 5,6, ...10) give that C = aD b 7 
with a = 0.27(5) and b = 0.92(13). 

The above derivation of the decay law (||) refers to the 
regime of quantum localization of chaos. For the inte- 
grable case one still has the rates T n ~ |c n | 2 ~ e - n / e eff 
where £ e f / is an effective length determined by tunneling 
in the classically forbidden region. However here the fluc- 
tuations in time are stronger since £ e ff depends on the 
local structure of invariant curves and islands in the intc- 
grable domain. Nevertheless the global decay at K = 0.5 
is in agreement with the 1/t law (see Fig. 2). In both 
cases (K = 7; 0.5), the quantum decay proceeds in a 
much slower way than the corresponding classical one. 
Due to the above reasons even in the case of mixed phase 
space the decay follows the 1/t law in accordance with 
numerical data of Fig.l. The 1/t behavior can continue 
up to a time t max oc exp (N/£ e ff) ~ l/r m j„ which is de- 
termined by the minimal decay rate in the system. For 
t > tmax the quantum survival probability P q (t) decays 
exponentially with time P q (t) oc exp (—t/t max ). 

The ionization rate T\ of an eigenstate ip\ localized at 
a distance n from the absorbing boundary is proportional 
to T\ oc exp(—n/£ e ff), therefore the number of such 
states is proportional to the measure /1 ~ n ~ In 1/r. As 
a result the probability to find a value T in the interval 
dT is dW/dT ~ d^,/dT ~ 1/r. Our numerical data for 
dW/ dT obtained in the localized regime confirm this 1 /Y 
dependence (see |l7j]). 

An interesting question is at what time scale the 
QPR start to deviate from their classical behavior. To 
test these deviations we determined them in two dif- 
ferent ways. The first one is defined as the time t q at 
which the ratio of classical over quantum probability is 
P(t q )/P q (t q ) — 0.9. This condition gives the first quan- 
tum deviation and was studied in [ jl6[] where it was found 
t q oc v^V oc y/l/h. This time scale can be explained on 
the basis of analysis of eigenvalue fluctuations of the evo- 
lution operator U in the complex plane or as the result 
of weak localization corrections |l6| . Our numerical data 
shown in Fig. 3 confirm this dependence even if the sit- 
uation is qualitatively different from that in |lq ], where 
the integrable component was absent, while here we have 
power law decay in a mixed phase space (see Fig.l). 

However at time t q quantum probability only starts to 
deviate from classical one while the behavior (||) can set 
in after a longer time tn- We determine this time in two 
ways: i) tn ~ In P(t q )/ In P q (t q ) = 0.9 (circles and stars 
in Fig. 3); ii) as the crossing point between the two extrap- 
olations for the classical behavior P(t) oc 1/t 2 and the 
quantum one P q (t) oc 1/t (crosses in Fig. 3). The data for 
stars in Fig. 3 were obtained by averaging over 50 values 
of N near a given value. We used two methods to deter- 
mine the time scale tn because the first definition gives 
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large oscillations related to the oscillations of the classical 
P(t) in lnt while the second definition allows to smooth 
out the effect of the oscillations. The numerical data are 
in agreement with the dependence tn ~ N cx 1/h. This 
time can be interpreted as the Heisenberg time scale, 
which is determined by inverse level spacings, and after 
which the quantum behavior becomes qualitatively dif- 
ferent from the classical one. 

The above scales t q and tn can be also seen in the 
quantum evolution of the Husimi function obtained from 
the Wigner function by smoothing over the size of fi cell. 
This evolution is presented in Fig. 4. For short times 
t < t q — 150 the classical and the quantum probability 
distributions in the phase space (n, 9) look rather sim- 
ilar. The quantum Husimi function reflects the under- 
line fractal structure of the classical distribution. For 
larger times the classical recurrencies are determined by 
more and more fine scales in the phase space and the 
classical distribution becomes localized around small is- 
lands near critical invariant curves. On such long times 
t > tjj = 3000 the quantum evolution in phase space in 
proximity of critical invariant curves becomes influenced 
by quantum interference effects that results in a qualita- 
tively different probability distribution (Fig. 4) . The pen- 
etration of quantum probability in phase space on smaller 
and smaller scales proceeds via a very slow tunneling pro- 
cess and as a result the Husimi function remains almost 
unchanged when the time is increased almost by two or- 
der of magnitudes (Fig. 4). 

At very long times t > t max the quantum decay be- 
comes exponential; this corresponds to eigenstates which 
are localized in the center of the island (Fig. 4). We have 
been able to observe such localized asymptotic Husimi 
function for N — 3 4 (Fig. 4). However, already for N = 3 5 
the time t max is too large to reach it in our numerical sim- 
ulations where we followed the evolution up to 3-10 . Ap- 
parently this time is already comparable with t max and 
this explains the strong oscillations of P q (t) for t > 10 6 
(Fig.l). At this time the quantum probability penetrates 
inside four small islands while at t » 1.5T0 7 the probabil- 
ity becomes concentrated inside the main central island. 
We would like to stress that t max grows exponentially 
with N (N ^ 1/K). As the result, the quantum decay 
1/t in semi-classical region proceeds up to enormously 
long times. In this regime the correlation functions prac- 
tically do not decay because C(t) cx tP(t) w const. 

In conclusion, we established a new universal law for 
probability and correlations decay in quantum systems 
both in the exponentially localized regime and in regimes 
in which classical dynamics has a mixed phase space. It 
should be possible to observe this universal behavior in 
the kicked rotator model which has been studied in ex- 
periments with cold atoms 



effect the fractal dimension of conductance fluctuations 
in soft billiards studied in |l0|]i"l| ]. According to these 
results the fractal dimension of these fluctuations should 
be cr = 2 — p/2 = 1.5. This regime should start at large 
times t > tji which requires to analyze conductance fluc- 
tuations with very fine resolution of magnetic field. The 
above value of a has not been seen in experiments |lC| ] 
that apparently indicates that the experimental resolu- 
tion was not yet sufficient to detect quantum Poincare 
recurrences at large times. 



13 1 . Similar effect should 



be also observable in the experiments with micrometer 
droplet lasing discussed recently in jl9| ]. The quantum 
decay (^) of P q (t) has exponent p — 1 that should also 
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FIG. 1. Classical (thick continuous line) and quantum 
(thin continuous and dotted lines) probability decays for 
K — 2.5 and N/k — 4. The two dashed straight lines show 
slope 2 and 1. The quantum curves correspond to N = 3 P 
with p increasing from 5 (upper continuous) to 10 (lower 
dotted). The starting conditions are two symmetric lines at 
n = ±iV/3, both for the classical and quantum evolution. 



FIG. 3. Dependence of quantum time scales t q (full dia- 
monds) and tH (circles, crosses, stars) on the system size N 
(see the text). The dotted and continuous straight lines show 
the theoretical slopes 1/2 and 1, for t q and tu respectively. 




FIG. 2. Classical (continuous lines) and quantum (dashed 
lines) probability decays for K = 0.5 (grey lines) and K — 7 
(black lines); parameter k = 5. The evolution starts at n — 
and the probability is absorbed for n < and n > 500. 
The two thin straight lines have slope one, while the dot- 
ted grey curve shows the fit for the classical decay at K = 7: 
P = 0.11 exp (-ft), 7 = 6.4 ■ 10 -4 . 
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FIG. 4. Left column: Husimi function at time 10 2 (top), 
5 ■ 10 3 (middle) and 3 ■ 10 5 (bottom), for the case in Fig.l with 
N = 3 8 . Right column: classical density plot of orbits surviv- 
ing up to time 10 2 (top) and 5-10 3 (middle). The right bottom 
shows Husimi function for N = 3 4 ,i = 10 7 . Distributions are 
averaged in a small time interval St = 20 near the given t 
values. The size of the phase region (8,n) is: < 9 < 2tt, 
—N/2 < n < N/2. The color is proportional to the density: 
blue for zero and bright red for maximal density (the color 
scale is the same for classical and quantum cases). 



